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Abstract 

Nonperturbative polaron variational methods are applied, within the so-called particle or worldline 
representation of relativistic field theory, to study scattering in the context of the scalar Wick - 
Cutkosky model. Important features of the variational calculation are that it is a controlled 
approximation scheme valid for arbitrary coupling strengths, the Green functions have all the cuts 
and poles expected for the exact result at any order in perturbation theory and that the variational 
parameters are simultaneously sensitive to the infrared as well as the ultraviolet behaviour of the 
theory. We generalize the previously used quadratic trial action by allowing more freedom for off- 
shell propagation without a change in the on-shell variational equations and evaluate the scattering 
amplitude at first order in the variational scheme. Particular attention is paid to the s-channel 
scattering near threshold because here non-perturbative effects can be large. We check the unitarity 
of a our numerical calculation and find it greatly improved compared to perturbation theory and 
to the zeroth order variational results. 


1. Introduction 


Variational techniques are widely used in the study of quantum mechanical systems. Most of these 
applications are based on the Rayleigh - Ritz variational principle for the energy of the system and 
physical intuition plays a crucial role in the construction of the variational wave function which is 
used. Applications of similar techniques to relativistic quantum field theories are, in contrast, rather 
limited. Intuitive arguments are hampered by having to deal with infinitely many degrees of freedom 
and the problem is further complicated by the need for renormalization. However, a well known 
example of a non-relativistic field theory where the variational approach has been highly successful 
is the polaron problem in condensed matter physics ii- Here it yields results within a few % 
compared to exact Monte Carlo results |^]. The main idea on which this success is based is the exact 
integration of the phonon degrees of freedom, thus obtaining an effective one-body system where an 
intuitive construction of a trial action is possible. The same approach has been applied to the Walecka 
model to study the one- @ as well as many- nucleon system. 

The generalization of this method to relativistic theories is possible and forms the basis of the 
work described in this paper. It should be stressed that this is a variational method which is not 
related to the usual variational approach in quantum field theory, the latter using trial wave func¬ 
tionals of the field in order to obtain an approximation to the effective action (see @ and references 
therein). Exact elimination of the light degrees of freedom (analogous to the phonons in the polaron 
problem) remains a central idea in this application. Furthermore, it is useful to employ the particle 
or worldline representation [2, 7 - 10] for the remaining heavy degrees of freedom of the theory. This 
formulation employs particle trajectories, parameterized in terms of a proper time, rather than local 
fields. Through its use the relativistic theory bears a startling similarity to the non-relativistic polaron 
problem and so it is not unreasonable to expect that some of the success enjoyed by the variational 
method so far will be repeated. 

In a recent series of papers [11-14] (henceforth referred to as I - IV) we developed the variational 
treatment of a simplified field theory of two scalar interacting particles. The Lagrangian of this 


theory is of the Wick-Cutkosky type [IT, IT] with a finite meson mass. As such, it has the same type 
of interaction as Lagrangians relevant to physical situations, for example any meson-nucleon theory 
or, in the limit of zero-mass. Quantum Electrodynamics (QED). On the other hand, it avoids the 
complications of a gauge theory as well as spin, chiral symmetry, etc. so it is fair to say that its 
primary purpose is to serve as a ‘playground’ on which to develop the necessary techniques of dealing 
with the strong coupling problem. Indeed, it is for this reason that, even after more than 40 years 
since Wick and Cutkosky first used the model to look for the relativistic bound state in the ladder 
approximation, it retains its popularity [17 - 23]. The generalization of the variational method to a 
realistic theory, namely QED, will be presented in a forthcoming publication [^] (for a preliminary 


account see p^ ). 

The Wick - Cutkosky Lagrangian in Euclidean space is given by 
^ 


( 1 . 1 ) 


where Mq is the bare mass of the heavy particle <I>, m is the mass of the light particle cf) and g 
is the coupling constant. As in (I) - (IV), we work in the quenched approximation of the theory 
and shall concentrate on those Green functions which involve one heavy particle interacting with an 
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arbitrary number of light particles. The generating functional for these Green functions is obtained 
after integration over the light degrees of freedom and is given by (see (I) and (II) for details) 

Z[j,x](x d/^exp (-^Mg ) / 'Dx{t)exp{-Ses[x{t)]-S 2 [x{t),j]} . (1.2) 

Jo V 2 / Jx(0)=0 


Here the effective action 


= fdt\x\t)-^^ r dh f dt2 / ^ 

Jo 2 2 Jo Jo J {^'x) Q + 


(1.3) 


is formulated in terms of the heavy particle trajectories parameterized by the proper time t. The 
meson source term 


S 2 [x{t)J] = -g I d^yj{y) I dt 


d f d*q 


0 J (27r)^ 


(1.4) 


produces n external mesons upon differentiation of Z[j,x] with respect to the meson source j{y) and 
setting j{y) = 0. 

The variational method due to Feynman is based on the construction of a suitable trial action St 
containing variational parameters and on Jensen’s inequality (or stationarity for complex actions) 




(1.5) 


where 

AS = Seff - St (1.6) 

and the averaging is done with the weight function e~^*. This averaging must be done exactly and 
therefore the trial action can depend at most quadratically on the function x{t). 

Several remarks are in order here. First of all, one may wonder in what way a quadratic trial 


action could in any way be a reasonable approximation to the effective action in Eq. (1.3) which, after 
all, is exponential in x{t). The important feature of Seg (which, by construction, is shared by the 
trial action St) is that it is non-local in the proper time t: For an ‘average’ path x{t) the separation 
x{ti) — x{t 2 ) tends to be correlated with the proper time difference ti — t 2 and so an unrestrained 
increase due to a quadratic term [x(fi) — x{t 2 ) in the trial action may be compensated for by a 
‘retardation function’ in ti — f 2 - Secondly, it is important to note that the approximation in Eq. o 
is a correctable one: if one so desires one can calculate arbitrarily precise corrections to the variational 
result by adding appropriate powers of AS on the right hand side of this equation. 

Another important feature of Eq. (|L^ ) is that it is easy to ensure that it reproduces the one loop 
perturbative result for any trial action whatsoever (see Ref. (I)). It goes without saying, however, 
that the more of the physics the trial action contains, the closer the numerical results of the ‘leading 


order’ result in Eq. (lA) should be to the true situation. In particular, in the present paper we shall 
examine meson-‘nucleon’ scattering near threshold (for simplicity, we refer to the heavier particle as 
the ‘nucleon’ and the lighter as the meson). One expects large distance (time) effects to be important 
at threshold since the particles move very slowly relative to each other, but, on the other hand, 
short distance effects also remain relevant because of the normal divergences of the relativistic held 
theory. In other words, particle scattering near threshold is one situation where it is important that 
the trial action contains information about both the infrared and the ultraviolet behaviour of the 
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theory (another application where this would be important is in a confining theory like Quantum 
Chromodynamics; see p^ ). In the variational approach under discussion it is once again the non¬ 
locality in the proper time which allows one to incorporate the physics of both these regimes into the 
trial action simultaneously through suitable behaviour of the retardation function for small and large 
values of — t 2 - 

It is the purpose of the present paper to examine to what extent this can be done for scattering, 
both analytically and numerically. A precursor to this work may be found in Ref. (III). There scatter¬ 
ing was discussed at zeroth order in the variational approximation (i.e. S’efr was just replaced by St)- 
Although a drastic approximation, resulting in a Gaussian dependence on all external momenta, the 
information contained in the zeroth order variational scattering amplitude was already considerable: 
for example, the forward scattering amplitude had an analytic structure which contained the cuts 
associated with the production of an arbitrary number of mesons in the intermediate state. In fact, 
this is another general feature of this variational method: the untruncated Green functions explicitly 
contain parts of all Feynman diagrams and their associated analytic structure of any order in the 
coupling constant (see (IV)). Furthermore, it is possible to perform the non-perturbative truncation 
of this Green function analytically, yielding a rather concise formula for Green functions involving an 
arbitrary number of external mesons. In that paper, the vertex function of the theory ( n = 1) was 
examined in detail and it was found that there are large non-perturbative effects, particularly near 
the pair production threshold. In the present paper we generalize the variational ansatz for the trial 
action further without modifying the variational equations for the on-shell particles. Thus the vertex 
function remains unchanged but the variational scattering amplitude ( n = 2 ) gets more freedom to 
approximate the true off-shell propagation. 

As in (III), we shall again pay particular attention to the question of unitarity near the lowest 
threshold of the forward scattering amplitude. Because unitarity is not antomatic in the formalism, 
it serves as a useful guide to its success. At zeroth order, unitarity is already satisfied to within a 
factor of two. We shall find that, as expected, this is improved in the present calculation (i.e. at ‘first 
order’ in the variational calculation). Guides such as this are particularly important for the theory 
under consideration as the Wick-Gntkosky model does not have a strong coupling limit |^, ]T| and 
so it is difficult to assess the quality of the variational calculation in the regime in which perturbation 
theory fails Q. 

This paper is organized in the following manner. In the next Section we describe the variational 
approach. In Section 3 we give the result for the general (2 -|- n)-point function. We use this result 
in Section 4 to derive the scattering amplitude for two external mesons (i.e. n = 2) as a general 
function of the Mandelstam variables s, t and u and examine its general structure as well as its region 
of convergence. In Section 5 we specialize to the s-channel near the first threshold and perform the 
necessary analytic continuation. Furthermore, we derive here important analytical results for the 
imaginary part of the scattering amplitude near threshold and in Section 6 we numerically investigate 
the question of unitarity. Finally, in Section 7 we conclude. 


^This should not be taken to mean that there are no interesting non-perturbative effects even at modest couplings - 
the qualitative and quantitative differences between the perturbative and variational form factors near threshold observed 
in (IV) are but one of these. 
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2. The Variational Approach 

The most general solvable trial action is of the form 


St[x]= [ dt^x^+[ dti [ dt 2 fih - t 2 ) [x(ti) - x{t 2 )f , (2.1) 

Jo 2 Jo Jo 

where f{ti — t 2 ) can be thought off as a ‘retardation function’ depending on the proper time lapse 
ti — t 2 between emission and re-absorption of mesons by the nucleon. 

In actual calculations it is preferable to work in Fourier space and write the trial action in the 
form 

OO 

St[b] = Y,Akbl ( 2 . 2 ) 

k=0 

where the bk are the Fourier components of the nucleon path 


, , t ^ 2^ , . 

x{t) = X - + > -— bk sm 


. / knt 

m - 

V P 


(2.3) 


Defining bo = x/v^2^, the free action is simply Sq = ^he nucleon propagator is given by 


1 rcc, 1 / /? r^\ 

^ exp --m 2 - — exp(- < A5 >s. 


(2.4) 


^ / Vb 
^ JVb e-^o 


It is well known that near the pole it is the Fourier transform of G 2 {x) which has a simple form. For this 
reason it is advantageous to apply the variational principle directly to G 2 {p) = f d^x exp(ip ■ x) G 2 (x) 
rather than to G 2 (x) using a trial action which is also a function of p ■ x. In (I) a simple trial action 
was constructed for the propagator by allowing an additional variational parameter A which can be 
thought off as a rescaling of the momentum p or a ‘velocity’ parameter. The form used was 


St = St — iXp ■ X 


(2.5) 


It was shown in (I) that in order to recover the correct pole behaviour of the propagator the proper time 
P must tend to infinity. In this limit all discrete sums over the Fourier modes Ak turn into integrals 
over the profile function A{E = k'xjP'). The ‘profile function’ A{E) is linked to the retardation 
function f{a) through 

e roo 

A{E) = 1 + ^J^ daf{a)sm^— . (2.6) 

The functional form of the retardation function, f(ti —t 2 ), may be left free and is best determined 
using the variational principle. In practice, however, it is more convenient to use a particular param¬ 
eterization containing a few variational parameters rather than a variational functional. This allows 
one to obtain analytic rather than numerical results. In our present numerical scheme it is also neces¬ 
sary for performing the analytic continuation from Euclidean time (where the variational calculation 
is done) to Minkowski time (where the scattering amplitude is evaluated). For the nucleon propagator 
it was found in (II) that it was quite sufficient to use trial actions which have a similar behaviour to 
the true action for small time intervals. For scattering, especially near particle production thresholds, 


4 



we also need the correct infrared behaviour of the theory. A particular example of a quantity which 
is sensitive to large distance scales is the imaginary part of a forward scattering amplitude near a 
threshold. It must rise linearly with momentum because of the limited phase space available to the 
total cross section which corresponds to it via the optical theorem. In Ref. (Ill) we used progressively 
more refined profile functions in order to investigate this. It was indeed found that whether or not 
the variational forward amplitude exhibited this linear rise (at zeroth order) was critically dependent 
on the infrared behaviour of the associated trial action. The variational principle came to our aid 
by providing the best form for the retardation function for all values of a and therefore we could 
determine what the small as well as the large a behaviour of /(cr) should be. In (III) it was shown 
that one can incorporate the correct large- and small-cr behaviour by using the following functional 
form for the retardation function (‘extended’ parameterization): 


Cl r - 


f(a) = -i - C^V^e 


—W2cr'\ 




with the requirement that wi > W 2 and C 2 > 0. The associated profile function is 

2Ci 


A{E) = 1 + 


E'^ 


E 

( E^\ 

E arctan- wi In 

1 H- 2 ) 

Wl 

V ^ 1 / 


(2.7) 


( 2 . 8 ) 


—4C'2\/^ 


(W 


W2 + E"^ + W2] - 


having branch points at E = iw 2 and E = iwi. The variational parameters Ci, C 2 , wi and W 2 were 
determined from the variational equations for the 2-point function G 2 {p) with the trial action of 
Eq. (|2.5D following the procedure described in (I) and (II). 


3. The (2 + n) point function 


The (2 -|- n)-point function can be derived in a similar way to the propagator G 2 {p)- In (IV) it was 
shown that the variational expression for the (2 -|- n) function is given by 


G'2,n(p,7'';{9}) oc / dl3 




W_9 I dr, 


. 2=1 


exp 




-<AS>a 

e 


d^x / Vxe-^d^d)] 


(3.1) 


In terms of the Fourier components the effective action in momentum space is 


S = 5[a;(t)] -|- i 



n+1 cxD 

bo • ~J 2 q 

i=l k=l 


(n) 

k 



(3.2) 


where 


<3l”’ = 


Tik 


E 

2=1 


Sin 


klTTi 

~T 


Qi 


(3.3) 


These additional terms on the right hand side (or, rather, their exponential) are just the plane waves 
appropriate for the n external mesons. The Euclidean momenta and the integration variables Oj are 
defined in Fig. ||. 0 

^ The unrestricted integration over all n in Eq. (1^ may be replaced by an integration over n as indicated in Fig. 
provided that all permutations of the meson momenta are summed over. In this way one wordline expression contains 
many different Feynman diagrams 
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Figure 1: The definition of the relative times and momenta. Note that tq = 0, Tn+i = f3, pi = —p 
and pn+i = p'. 


Given a true action of the above form, it was a natural choice in (IV) to take the trial action for 
the n-point function to be of the form 


St = St + i 


n+l 




2=1 


k=l 


(3.4) 


where A and are variational parameters which turned out to have the same variational equations 
independently of the number of external mesons n. ^ 

Actually, one can do somewhat better than this without much effort and, indeed, it turns out to 
be physically sensible to choose a somewhat more general trial action than Eq. (p.4|). The reason is 
the following: although we only consider on-shell Green functions, the intermediate nucleon lines in 
Fig. g are in general still off-shell. Therefore it might be useful to associate a different A with each 
of the (n -|- 1) nucleon momenta in this action, thus allowing the trial action additional freedom to 
adjust itself to the fact that these internal lines have different virtualities. As we shall see, for the 
correct position of particle production thresholds in the intermediate state this choice of trial action 
will prove to be useful. 

In the following, we briefly summarize the relevant results for a trial action of this form. The 
discussion follows (IV) closely, where the interested reader may find further details. Grucially, as we 
shall see shortly, the variational equation for the A’s associated with the external legs reduce to the 
same equation as one found for the nucleon propagator, resulting in the same mass renormalization 
as before. Therefore, a consistent truncation of the Green function remains possible. 

In short, the trial action which we shall use is given by 


St = St + i 


n+l 


bo ■ X - X Qk"'’ ■ bk 


2=1 


k=l 


(3.5) 


^The reason why the combination \Au appears in the last term of Eq. (3.4) is connected to the truncation of the 
Green function - see Ref. (IV). 
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where 


/^(n) 2y^ . 


Z=1 


klTTi 

~T 


{^i+lPi+1 ^iPi) ) 


(3.6) 


and St is given by Eq. (^). Note that in the limit in which the coupling goes to zero the true action 
S becomes 


So — So + i ^ 

i=l k=l 

so that Ak and Aj must go to 1 in this limit. Writing the true action as 


n+l 


-bo ■ am - • bk 


(3.7) 


S = So + Si (3.8) 

one obtains for the path integral of the exponential of the trial action, up to irrelevant overall constants, 


• - 00-1 
Vxe-^^= 

.k=0 


exp 


1 


R °° 


(n)= 

k 


k=l 


where 


Vx = / (Tx 


rx{0)=x 

lx{0)=0 


Vx . 


Further, as part of < AS" > we need 


oo Qinr 


<So-St >s^- — -^(l + ^o) + 2^ - l) ^^^(1 + Afc) + Y 


oo Q(n) . qH 


(3.9) 


(3.10) 


, (3.11) 


k=l 


k=l 


For brevity we have defined the quantity P to be sum over all the nucleon lines’ momenta weighted 
by the respective proper time intervals, i.e. 


n+l 

p = Y ■ 

i=l 

P is the equivalent quantity containing the relevant factors of Aj, namely 

n+l 

P =Y ■ 

i=l 


Combining these terms we obtain 

— < §0 — St >a / Vxe~^* = exp 


exp 


1 


n+l 


2 E 


2) 


k=0 ^ ' i=l 

Finally, the weighted average of the interacting part of the action is given by 

-.2 rS 1 /-I / 


< + >+= - 


g 


IGvr^ 

where e(x, y, z) is defined to be the exponential 


/ dtidt 2 [ dxi e(mjl{a,T), 

Jo ppa,T) Jo V 


il{a, T) 


, Xi 


fi{x,y,z) = exp 


x"^ 1 — z \ 

-z 

2 z 2 ) 


(3.12) 


(3.13) 


. (3.14) 


(3.15) 


(3.16) 
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The quantity /U^(cr, T) appearing in the above equations is defined in terms of the ‘pseudotime’ T), 


i.e. 




a 


AqP 

AoP 




+ 


8/3 


TT^ 


E 


sinE^'l . 


V /3 ;V 2 / 3 ; 


(3.17) 


Here we have introduced the relative and total times a = ti—t 2 and T = {ti+t2)/2 and the pseudotime 
/r^, after taking the /3 ^ oo limit, becomes 

4 


1 sm‘^{Ea/2) 

dE 


" ; l ME) — & — ■ 

The analytic structure as well as the asymptotic behaviour of the pseudotime in Eq. ( 3.18| ) is of great 
importance. This analytic structure is closely linked to the analytic structure of the profile function. 
In particular, parameterizations of the profile function which exhibit the same ultraviolet behaviour 
as the true action lead to a pseudotime which grows linearly with a for both c —> 0 and c —> oo 

and does not have poles or cuts in the half-plane where TZe (a) is positive (see III). 

The quantity IT^”^ is defined by 


(T ^ 2^^ 1 


IEW = -p - 

(3 TT 


oo 

TeS 

k=l 


kirT 


k 


cos 




/ kTra 


sm 


V 2/3 


It is equal to 


1 


IHW = _(Ai Pi + A„+i pn+i) + {\Ti - ti\ - In - t2\ ) (Ai+i pi+i - Xi Pi) 

i=l 


which simplifies to 

IT^”^ = sign(cj)| 


Tfe+i - min(ti,t2) 


Afe-l-l P6-I-1 


(3.19) 


(3.20) 


(3.21) 


a—1 

+ E/ + 

i=b-\-2 


max(ti,t 2 ) - Ta-l 


XaPa } 


= CrXaPa 


a > b + I 
0 = 6+1 


For the case a = 6 + 2, the sum appearing in Eq. ( 3.21 ) is defined to be empty. In other words, 
is equal (up to a sign, which is not relevant as only the square of IT^”) enters) to the integral of the 
proper time multiplied by the modified nucleon’s four-momentum Xp for the duration of the exchange 
of the internal meson. Particular cases of interest, for the purpose of the eventual truncation of the 
Green function, are those where ti ^2 are both smaller or both larger than all p, in which case = 

aXipi and aX^+iPn+i, respectively. 

Having calculated the above averages, one can now write down the expression for the (untruncated) 
(2 + n)-point function to first order in the variational framework 


No 


' n+1 


G'2,n(p,p';{g}) = ^ E ills' 

S' I i=i 


Svr^ 


dai exp 
r*i 


(Xi 


~~2 “ (1 “ 


(3.22) 


/ dti [ dt 2 ' / dxi e\mp{a), 
Jo Jo P \^) Jo \ 


-il+W 

p{a) 


Xl 


■ exp 


























Here Nq is defined to be 


1/1 
A'^o = -TTTTT exp 1 — 


AiO)J - 

For i = 1 and i = n + 1 we get the same on-shell relation between the bare mass and the renormalized 
mass as before, e.g. for i = 1 


M« + 2 + I? /f C 


Xl 


(3.24) 


where the quantity H has been defined in Ref. (I). Further, we get the same variational equations for 
these two lambda’s as before and so we write 


Ai — A„+i — A 


(3.25) 


The truncated (2 + n)-point function G^j^{p,p'‘,{q}), valid for on-shell external nucleon momenta, 
now reads 


G 2 ^n{P,P'-,{Q}) = If I / dai exp 


(Xi 


(M2(l-(l-A)2)+p2(l_(l_A,)2)) 


da ^ . XaM 

/ e m/r a), —^ , n 

/o ^ V p{(y) 


dtidt2 ( , , ^ 

e mp{a), -— , u 


//2(cj) p{a) 

where Viqi} denotes all possible permutations of g*. Eq. (|3.26|) is the main result of this section. 


(3.26) 


4. The (2+2) scattering amplitude 

In this section we focus our attention to the description of the scattering amplitude for meson-nucleon 
scattering. From the last equation of the previous section, substituting n = 2, we immediately obtain 

G2,2{p,qi,q2) = ff- [ da2 expi-^\M^{l-{1-\f)+ pI{ 1-{1-\2f) 


+ ^ '^{p,qi,q2]a2) 


(4.1) 


The function H is defined as 


'■‘{p,qi,q2;a2) = / dxi 

Jo 


dtidt2 ( , , '' 

e mp[a), — ^ , xi 


+ (cj) 


p{a) 


f°° da f i \ XaM 

/ -wrw(^‘2e[mp{a),——,xi 

lo +(o-) \ p{a) 


(4.2) 


The ti — t 2 integration region in Eq. (T2) is shown in Eig. Depending on the ordering of the 
internal meson’s emission and absorption times t 2 and ti with respect to the proper times at which 
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the external mesons interact with the nucleon, i.e. ri and T 2 , the function assumes a different 
value. Explicitly, 

^ t2 <Ti ,Ti<ti < T2 

= itl-t2)X2P2 Ti<ti,t2<T2 

= (t2 - t2)X2 P2 + {h - T2)X P3 Ti < t2 < T2 ,T2 < h 

= {Tl-t2)Xpi+a2X2P2 + {ti-T2)Xp3 t2 < Ti ,T2<ti. (4.3) 



Figure 2: The (ti - t 2 ) - integration region is shown enclosed by the solid line. The (dotted) triangular 
regions are those relevant to the dressing of the external legs and have been removed in the truncated 
amplitude. Also shown are the relevant diagrammatic representations (for the direct diagram) of the 
action in the particle representation 


Although the integration variable 02 in the expression for the Green function in Eq. ( 3.26 ) is a 
time difference, the integration variables ti and t 2 in Eq. (^) are still absolute times. This leads 
to a spurious dependence of and the integration limits of the (ti - t 2 ) - integration on the 

absolute times ti and T 2 ^ In order to eliminate this spurious dependence it is convenient to make 
a transformation of the integration variables ti and t 2 to two variables which are expressed as time 
differences. Clearly cr = ti — t 2 is a suitable choice for one of these. It is convenient to define the 
second variable X 2 in a different way for each of the integration regions depicted in Fig. We shall 
discuss each of these regions below. 

Before proceeding, let us make the change to the Minkowski space Mandelstam variables s, t and u 
(as well as Minkowski and q'^ for the external mesons), which are defined in terms of the Euclidean 


should be remembered that n and as have already been integrated over in Eq. (3.26) in order to truncate the 
Green function, so a dependence on these variables in Eq. (3.26) would be non-sensical. 
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momenta pi, ps, qi and q 2 through 


s = -{pi+qif 

t = -{Pi-Psf (4.4) 

u = -{Pi + q2f ■ 

In contrast to the usual formalism this change from Euclidean to Minkowski four-momenta is a trivial 
one in the wordline approach, as there are no four-momentum integrals left to be done in the expression 
for the Green function. 

It is rather useful to break up the function H not only into the components corresponding to the 
two vertex corrections, the propagator correction and the box diagram depicted in Fig. i.e. 

E;(s, t, q^, a 2 ) = '^v{s, q^] 02) + Sy (s, a 2 ) + Hp(s; 02) t, g^, g'^; 02) , ( 4 . 5 ) 


but also to immediately extract their leading 02 - behaviour. The reason for this is that Q !2 is the 
proper time conjugate to the intermediate nucleon’s momentum and so one might intuitively expect 
that the large a 2 - behaviour is particularly important close to the threshold, i.e. close to where real 
particles could be produced in the intermediate state. This is analogous to what happens for the 
nucleon propagator itself - as indicated in Section 2, the on-shell propagator corresponds to letting 
/3 —> 00 . This physical interpretation associated with the use of the Fock-Schwinger proper time 
formalism makes it advantageous over the use of the usual Feynman parameters. 

Explicitly, the component of S corresponding to the bottom left hand integration region in Fig. ^ 
may be written as 

“;y(s,g^;a2) = ^^(s,^^)+ Sy(s,g^;a2) , (4.6) 

where Fy(s, g^) does not depend on q ;2 and Fy(s, g^; 02 ) goes to zero as Q !2 —> 00 . In these variables, 
the constant piece Hy(s, g^) for the direct diagram (we shall discuss the crossed diagram later) becomes 

Wv{x2),x^ . (4.7) 


^^(■5,9^)= [ dxidx2 
JO 


da 


a 


p^ia) 


mp{a] 


a 




Here the variable X 2 is defined as the time interval ri —12 between the external and internal meson’s 
interaction time with the nucleon, scaled by a, and Wy{x 2 ) is given by 


Wy{x2) = X2A2 + X2^(A - A2) AM^ A 2(1 - X2) s (A2 + X2(A - A2)) - Ax2g 


(4.8) 


On the other hand, the piece of Sy which depends on a 2 is 

Ev{s,q‘^-,a 2 ) = - [ dxidx 2 [ da ^ ^ e (mp{a + « 2 ), Wv {—— X 2 ),xi 

Jo Jo p^{a + a 2 ) V p{a + a 2 ) a + a 2 

(4.9) 

where a has been shifted by 02 and X 2 has been rescaled accordingly. Also, it is easily seen that the 
expression for the second vertex (top right hand region in Fig. is the same, except that g^ needs to 
be replaced by g'^. 

In a similar way to the vertex function discussed in (IV) the representations Eqs. (4.7) and ( ^.91 ) 
for Hy (s, g^; 02 ) and Fy(s, g'^; a 2 ) do not necessarily converge for all s. In fact, if we restrict ourselves 
to the case where the external mesons are on-shell (i.e. g^ = q'^ = mJ) then Eqs. (^) and (O) are 
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only valid representations if s > 0 . | If one is interested in s < 0, then Eqs. (^) and (|4.9|) must be 
analytically continued through the use of a contour integration in the complex cr - plane. 

It is worthwhile at this point to remind oneself of the corresponding situation for the vertex 
function (see Ref. (IV)). Here the original (real) representation for G 2 ,i(g^) was valid for all values of 
below the physical region for nucleon pair production (i.e. < 4M^). For larger than this the 

same contour rotation as above had to be performed, resulting in a cut and hence an imaginary part 
for the amplitude, starting at the nucleon pair production threshold. The situation for Hy is reversed, 
as we have just indicated that this is real for s > 0 (and thus real in the s-channel - see Fig. |3|), while 
the corresponding quantity seems to become complex in the t and (most of the) u channel. At first 
sight this seems fatal ~ after all, the contribution from the direct diagram should have a cut beginning 
at s = (M + m)^ and should be purely real in the u channel. We shall address this issue and its 
solution below, after having discussed the contributions to S from the other two integration regions. 



Figure 3; The Mandelstam plot for the scattering process 

The part of H originating from the propagator correction may be decomposed in a similar manner 
to 'Ey. This time, however, there is also a term which goes linearly with 0 : 2 : 

^p{s\ 02 ) = —Oi2 ‘=‘p('5) + '^p('S) + 0 : 2 ) • (4-10) 

There is a good reason for the existence of this linear term for this particular integration region; if we 
would expand our expressions in a perturbation series (so that the diagrams in Fig. ^ would actually 
correspond to Feynman diagrams), this term would be responsible for the self-energy correction S(s) 
of the intermediate off-shell propagator. Mathematically, the term linear in a 2 results from the fact 
that the integration region in the complementary variable to a (which may be performed analytically) 

®This requirement originates in the fact that grows like a for large cr, so that for the cr - integral to converge 

we require Wy > 0 for all 0 < 2:2 < 1. 
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grows like 02 , as may be seen in Fig. Clearly this is the only region for which this happens. It is 
also no accident that the expression for Hp looks very similar to those encountered in the discussion 
of the renormalization of the nucleon propagator in (I) 


^P 



da 

T^) 


a 


e ( m^ia), AM, xi 
fi{a) 


e ( mix{a) 


a 




■A2\/i,xi 


(4.11) 


the only difference being that there it only applied to the on-shell propagator, while Hp also contains 
off-shell information. Indeed, as was proven for the general n-point function in (IV), the variational 
calculation - if expanded in the coupling - reproduces perturbation theory to one loop order exactly, 
including the off-shell behaviour of the internal propagator. 

The constant piece contributing to Hp is given by 



(^mn{a), 


a 

fj.{a) 


X2^/s,Xl 


(4.12) 


while the piece that vanishes as 02 00 becomes, after shifting the a integration by a 2 as before. 


zp[s-,a2) = / 
Jo 


(4.13) 


Clearly, the region of convergence of the representation for Hp (Eqs. ( |4.10[ - |4T^ )) is the same as that 
of Hy, i.e. s > 0 . Also, it is relevant to the subsequent discussion of the region of convergence of the 
02 - integration that in the s channel Hp is positive. 

Finally, the integration over the region where the action is represented by the box diagram is given 

by 


EB{s,t,q^,q'‘^-a 2 ) = [ dxidx 2 [ da ^ - -e (m^j.{a + 02 ), ^ . IFp, xi) , (4.14) 

Jo Jo ^i^{a + a2) \ ij{a + a2) ) 


where 


2 _ (o'A-|-a 2 A 2 ) . -v 71^2 


1F^ = 


(cj -I- ^2)^ 


a\ 


a ^ 02 


ao2XX2 


{aXM + 02 X 2 s)-[ - {1-X2)x2t--. —;-^ X 2 q + {I - X 2 )q . (4.15) 


(cr-|-a2)^ 


This time there is not even a constant term in 0 : 2 , and the representation converges (for on-shell 
external mesons) as long as t < 4M^, i.e. in the u and s channel. 

In short, the contribution to the scattering amplitude due to the direct diagram is given by 


02 


A“(s,t,u,g^g'^)=.F(s,g^g'^) da 2 exp <( ^^[1 - (1 - A 2)1 [5 - 5(s)] + ^H(s, t, a 2 )| , 

(4.16) 

(4.17) 


where 


s(s) = 


1-(1-A2)' 


(l-(l-A)W + £2 =p(») 


and 


T-l 2 l2\ ^ 

JF{s,q ,q ) = —— exp 


^ {^vis,q‘^)+^vis,q'‘^)+E%{s)^ 


(4.18) 


Furthermore Ni = XZjN^ where Z is the residue at the pole of the propagator. Note that for 
reasonably small couplings s is ‘close’ to M^ and that for s > M^ it is larger than M^. 
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Let us now come back to the question of the analytic structure of this amplitude. It is clear from 
Eq. ( 4.16| ) that not only does one need to address the convergence criteria of the a - integration, 
but also simultaneously those of the 02 integration. It is for this reason that the analytic structure 
of the scattering amplitude is much richer than that of the vertex function in (IV). It is clear from 
Eq. ( I.IC ) that for the 02 - integration to converge we require s < s (note that A is always less than 
1). However, s, using the representation for Hp in Eq. ( 4.11 ), is only defined for s > 0 . So at the 
moment our Green function is only defined in a region 0 < s < s and t < 4M^ , i.e. except for 
a limited region of the u - channel (Fig. it is only defined outside the three physical regions for 
the process. Importantly, however, for the analytic continuation there does exist a region where it is 
defined. Because of the complicated interplay between the convergence of the a - and Q !2 - integrations 
one needs to analytically continue the above representation on a case by case basis. In the present 
paper we are interested in doing this for the s-channel. 

Before we proceed, let us write down the amplitude for the crossed channel. It is given by 


A crossed 

^ 2,2 




A direct 

^ 2,2 


{u,t,s,q''^,q^) . 


(4.19) 


Hence, again for on-shell external mesons, the representations above are only appropriate for 
in the region 0 < rt < s and t < 4M^ . As may be seen in Fig. this includes some of the physical 
region of the s channel. Indeed, as we shall concentrate our attention on the region near the threshold 
at s = (M + m)^, this representation of suffices for our purposes. 


5. Unitarity in the s channel near threshold 


In order to be able to describe scattering of on-shell mesons off the nucleon in the s channel, the 
convergence requirement of the 0.2 - integration of the direct diagram requires the amplitude in 
Eq. ( |4.16 ) to be analytically continued so that it converges for s > s. This may be done in an 
analogous manner to a Wick rotation as in Refs. (HI,IV): for s infinitesimally above the cut, i.e. 
s s + ie, one may deform the - integration in such a way that it runs along the positive 
imaginary axis. This is possible because the contribution of the quarter circle at infinity vanishes and 
because the analytic behaviour of the pseudotime /^^(cr) is known. Changing variables to 02 ^ 0 : 2 /*) 
and suppressing the dependence on = vr? for notational simplicity, one obtains 


A direct 

^ 2,2 


{s,t) = 


iT{s) (iQ! 2 exp I - (1 --^ 2 )^] [s - s(s) + *e] + ^ 2(s,t;ia2) I (5.1) 


exp 




- 1 


where the latter expression, obtained by adding and subtraction a ‘modified Born term’, does not 
require an ie prescription for convergence and is therefore more suitable for numerical evaluation. 

This expression for the amplitude converges in the region of interest. Together with the crossed 
diagram, which does not have to be analytically continued as long as one does not move more than 
about one nucleon mass away from the lowest threshold (and therefore remains real), it yields the 
complete scattering amplitude. Furthermore, it should be noted that for the special case of forward 
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scattering the expressions for S naturally simplify somewhat. This is particularly true for the contri¬ 
bution of the box diagram to the action, where the X 2 - integration may be carried out trivially when 
t = 0 and = q'‘^. Also, the two vertex contributions are, of course, identical whenever the latter of 
these conditions is fulfilled. 

It is relatively straightforward to evaluate the crossed diagram numerically. Both the a 2 - as well as 
the a - integral in E{u,t; 02 ) converge exponentially, so simple Gaussian quadrature is quite adequate 
(we typically used 64, 72 or 96 Gaussian points per integration). The direct diagram, however is much 
more difficult to handle as in this case the integrands are oscillatory functions. In particular, the 02 
- integration only converges relatively slowly as 02 ^ 00 . Moreover, the separate contributions to 
S(s,t;a 2 ) arising from the vertex, propagator and box diagrams oscillate with different frequencies 
for asymptotic 0 . 2 , so that the numerical behaviour of the exponential of S(s,t;Q; 2 ) is rather erratic. 
This means that for a reasonably reliable result one needs to integrate up to rather large values of 
OL 2 - We do this integral using the adaptive routine DOlASF for sine and cosine transforms from the 
NAG Fortran library. 

We shall now turn to the question of the unitarity of the scattering amplitude near the lowest 
threshold. As mentioned in the Introduction, this is an interesting quantity to examine because one 
expects non-perturbative effects to be particularly strong at thresholds. Below the first inelastic 
threshold, unitarity implies that the total elastic cross section is related to the forward scattering 
amplitude. In the centre of mass system 

g'ei('S) = r > (5-2) 

/|p| VS 

where p is the centre of mass three-momentum and the elastic cross section is given by 

o-ei(s) = J dQ -^^^\A{s,t)\^ . ( 5 . 3 ) 

Near threshold the elastic scattering cross section is, of course, dominated by the real part of the 
amplitude because the imaginary part has to rise linearly with momentum from zero at threshold. 
This means that Uei is not particularly sensitive to the exact position of the threshold. For the 
imaginary part of the forward scattering amplitude, however, the threshold position is of supreme 
importance. As was seen at zeroth order in (III), a linear rise with |p| as well as a correct threshold 
position was not assured for an arbitrary trial action. In particular, the large-ci behaviour of the 
retardation function entering this action was found to be crucial. 

At first order one does not expect this to be quite as severe, as at least agreement with one- 
loop perturbation theory (which has the correct threshold behaviour) is assured for an arbitrary trial 
action. Nevertheless, the basic point remains the same: threshold behaviour is governed by the large 
distance behaviour of the action and therefore we expect some sensitivity to the large - a behaviour 
of the retardation function. Hence it is advisable to spend some time examining the dependence of 
the imaginary part on the form of the profile function A[E). 

5.1 Xmyl(s, 0) near Threshold 

Before doing this, it is actually very instructive to have a look at ImA(s, 0) in perturbation theory. 
To leading order in |p|, for each of the diagrams which includes a vertex correction, one obtains near 
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the threshold 


Im 0) = —32a^7r 




(4M2 — vn?‘){M + m)m‘^ 
the term involving the self energy of the intermediate propagator is 

T\/fi 

Im 0) = 32a^7r 


(2M + mY{M + 


and the box diagram yields 

Im 0) = 32a^7r 




P\ +0 (p2) 


+ O (p') 


+ O (p') 


(2M — mY{M + m)m? 

Altogether, the leading contribution to the imaginary part is therefore given by 


Im^pert.(^^O) = 128a27r 


(4M2 - m2)2(M + m) 


IpI + o (p') 


(5.4) 


(5.5) 


(5.6) 


(5.7) 


Even though each diagram separately is divergent as the meson mass m goes to zero, this divergence 
cancels in the complete amplitude. For the physical value of ^ ~ ^ means that the individual 
diagrams are something like 50 times as large as the final result. This makes the task of calculating 
the imaginary part of the variational expression for the forward scattering amplitude in Eq. G) even 
more formidable than it already is - not only does one have the numerical problems mentioned above, 
one also should expect large cancelations between the various components of H(s,t;a 2 ) ■ 

At first sight it may appear that the lowest threshold will open up as soon as s is large enough 
so that one needs to work with the complex representation in Eq. (|5.1| ) rather than the purely real 
representation of Eq. ( 4.16 ), i.e. s > s(s). Certainly, for 0 < s < s{s) the amplitude must be real, 
because in this region Eq. ( 4.161) is valid. It is useful to examine the reason why Eq. (|5.lD , which is 
also valid in this region of s, does not develop an imaginary part here. To this end we shall define 


and 


pW = 1[i-(i-A2)2] [s-s(«)| 


g 


G{a 2 ) = exp I — S(s, t; ^ 2 ) | - 1 . 

Hence the imaginary part of the forward scattering amplitude is related to 

poo 

/ da2 [cos{a2'D)TZe G{a2) — sm{a2'D)ZmG{a2)] 


(5.8) 


(5.9) 


(5.10) 


Using the facts that G(a 2 ) = G*(—a 2 ), that G(a 2 ) has a cut along the positive imaginary axis and 
that it vanishes sufficiently fast as 02 ^ 00 , one may easily prove, using dispersion methods, that 

poo poo 

/ da2 cos{a2'D) TZe G{a2) = — / da2 sin(Q:2|II|) Xm G(a2) • ( 5 - 11 ) 

40 40 

Hence the two contributions to the imaginary part of the scattering amplitude in Eq. ( 5.10| ) cancel 
each other when D < 0, i.e. s < s(s), as required. For s > s(s) they remain identical in magnitude, 
but add. This is a useful property as it means that the amount of computer time required in the 
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numerical integrations involved in the calculation of the imaginary part of Eq. (5.1) is cut by a factor 
of 2. 

Just because Eq. (ED may develop an imaginary part for s above s does not mean that it will do 
so. In fact, one may use a similar argument to that used in the investigation of the threshold position 
in (III) to show that the imaginary part only starts somewhat above s. Let us write H(s,0;a2) in 
Eq. ( 4.16| ) as a Laplace transform: 




Svr^ 


POC 

H(s,0;a2)= / dEp{E,s) 

JO 


,-Ea2 


(5.12) 


Here p{E, s) is essentially the counterpart of the weight function p{E) encountered in the zeroth order 
calculation in (III). It was found, in fact, that p{E) is only non-zero for E larger than a critical value 
Eq and this will be also true for p, as will be shown below. Expanding the exponential in Eq. ( 4.16 ) 
with s —> s -|- ie and performing the 0:2 - integration then yields 

TmH“(s,0) = 7r.F(s) ^ - / dEi...dEn ^Ei, s)...p{En, s) 5 (Y, Ei - V{s)] (5.13) 

Only when the J-function is fulfilled an imaginary part develops. With this observation we proceed 
to obtain an analytic expression for the leading behavior of the imaginary part of the amplitude near 
threshold. Given the difficulties with the numerical evaluation of Eq. (5.1), in particular in view of the 
large cancelations expected between Ey, Ep and Ep, the possibility of having an analytic expression 
is indeed fortunate. The key to obtaining this threshold behaviour is the observation that below the 
first inelastic threshold only the n = 1 term in Eq. ( ^.13 ) contributes. This amounts to expanding 
the exponential in the expression for the amplitude Eq. ( 4.16 ) and keeping only the first term; i.e. 

Im 0) = ^ E{s) Im / da 2 H(s, 0; 02 ) Eq < V{s) < 2 Eq . (5.14) 

’ ovr^ JO 

Below the first inelastic threshold the imaginary part of the forward amplitude is then simply given 
by 

TmAfY\s, 0) = .F(s) p ( V{s), s ) 0 {V{s) - Eq) (5.15) 

where we made use of the Laplace transform of E defined in Eq. ( ^.12 ). 

In practice, the weight function p{E, s) can be determined neither numerically (as inverse Laplace 
transforms are notoriously unstable) nor analytically because of the complicated 0.2 - dependence of 
the function E. The only case where the inverse Laplace transformation can be performed is near 
threshold which is known to be dominated by large values of 02 . In this case we can approximate the 
pseudotime by its asymptotic behaviour 

0 L 2 


p ( 02 ) 


a2^oc 


H(0) 


+ 4^(0) + O (e-"2) 


(5.16) 


This allows us to evaluate the corresponding asymptotic expressions for the E-functions. For example, 
the propagator diagram’s contribution becomes 


’(^ 2 ) 


H(0) 


da / dxi 


a 


exp 


a 2 + a + 4^(0) 

•exp —4 j1(0)^(0) 


1 — xi H(0) .0 \ , X 

H-—A2'S3 :i I ( 0:2 E a) 


\^2H(0) xi 

m? 1 — xi 

,2H(0) XI 


2 

H(0) 


A2 s Xl 


(5.17) 
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The above form is not a fully consistent expansion in inverse powers of 1/02 since the exponential 
function will also contribute to such terms. Nevertheless it is a convenient form to perform the inverse 
Laplace transform as we may use (Ref. Eq. (5.5.9)) 


C 


-1 


a2 + b 


■a 2 {E) = Q{E-a) 


to obtain 


Here 


pp{E,s) ~ 


^ ^ ^ (ixi 0 ( E — Eo(s, xi) ) exp —Q^(s, xi, E) ^(0) 

1 — xi H(0) 2 


87r2 E2 Jo 

Eo{s,xi) = 


2H(0) XI +^^2sa;i 


determines the minimal value of E (see below) and 

1 — xi H(0) 2 


Q"(s,xi,E) =4yl(0) 


2H(0) xi 


A 2 SX 1 + E — Eo{s, xi) 


(5.18) 

(5.19) 

(5.20) 

(5.21) 


is a kind of effective momentum transfer (one should remember from (III) that 6 ^(0) is the mean 
square radius of the dressed particle in lowest variational order ). 

Similarly we obtain for the asymptotic form of H in the box diagram 


^ 5 ( 02 ) ^ ^(0) da dxi 


a 


0 02 + 0 -+ 4^(0) 


exp 


-Eo{s,xi)a2 - Q‘^{s,xi,Eo) ^{0) 


• exp 


m? 1 — xi H(0 )xiA2 


+ 


^A(s + — m^) — A 2 S^ a | (5.22) 


2H(0) XI ' 2 

where again we have neglected terms of order l/a 2 iu the argument of the exponential function. Using 
the formula (5.18) we get 


Pb{E,s) ~ 


^2 .1 


87r2 


dxi Q { E — Eo{s, xi) ) 


Al(0) 


Ef{s,xi,E) 


exp 


-Q\s,xi,E)m 


(5.23) 


where 


Ei{s,xi,E) = ^ -—— + (^X{s + M‘^ - m^) - \ 2 S^ + E - Eo{s,xi) . (5.24) 


2H(0) xi ' 2 

Finally the H-function in the vertex diagram has the asymptotic limit 


Sy(Q!2) —> —^(0) / da dxidx2 


a 


a 2 + a + 4^(0) 


exp 


-Eo{s, xi)a 2 - Q^(s, xi, Eo) ^(0) 


• exp 


m? 1 — xi H(0 )xiA2 


+ 


(^Ax 2 (s + — m^) + (1 — 2 x2)A2s) a | .(5. 


2H(0) xi ' 2 

Consequently the corresponding weight function for the two identical vertex contributions becomes 


25) 


pv{E,s) ~ -2-^A{0) dxi dx 2 0 ( E - Eo(s,xi) ) exp [-Q^(s, xi, E) ^(0) 
87r2 Jo I 


m? 1 — xi ^ H( 0 )xiA 2 


2H(0) xi 


(^Ax2(s + - mJ ) + (1 - 2x2)A2s)^ 


-1 


.(5.26) 
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The X 2 -integral is the same as the one for combining denominators by the Feynman parameterization 


fo [0 + 6x2 


1 

ab 


and thus we obtain 

me , s) dx, e ( E - E„(s.xO ) 

Combining all terms gives the simple formula 

HE. X) = £^A{0) e { E - EM X.) ) (i - 


exp 


(5.27) 


. (5.28) 


exp 


-Q2(s,;ri,F;)e(0) (5.29) 


which is evidently positive. This equation, together with Eq. ( ^.15 ), directly determines the imaginary 
part of the forward scattering amplitude. Hence one obtains a positive total cross section, as required. 
In the following, we shall first examine for which s this imaginary part first appears and subsequently 
shall derive an analytical expression for the total cross section at threshold. 

5.2 Threshold Position 


Let us now determine the threshold position. Eq. (5.29) should hold in the vicinity of the elastic 
threshold. Approaching it, the range of the xi-integration shrinks due to the step function and is 
restricted to the interval xj^ ^ < xi < x^"*"^ where 


xl"^^ = 


1 


A(0) Ai 


E + 


m 


2x4(0) 


9 \ ^ 


(5.30) 


Exactly at threshold x^"*"^ and x^ ^ coincide and only the minimum point 

Xl = 


m 


A(0)A2\/s 


is left over. Eq. (5.31) determines the lowest threshold position to be 

.^thr. ~ 'P('5thr.) ~ -®0(®1)'^thr.) 


(5.31) 


(5.32) 


or 


Eq = Xmy/s — 


m 


2A{0) 


(5.33) 


which is the statement that that p{E,s) is only non-zero for E larger than Eq. Inserting Eq. (5.33) 
into Eq. (5.12), one finds the threshold condition 


n 2 A 2 - 


\ 

m 


1 -( 1 -A 2 )^ [sthr. - s(sthr.) ] 


(5.34) 


The parameter A 2 is undetermined up to now. Certainly, one might try to fix it via the variational 
principle, which, however, in practice would seem to be rather difficult. An alternative, which is what 
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we shall do here, is to fix it by demanding that the lowest threshold (i.e. n = 1) coincides exactly 
with A/ithj. = M + m. ^ 


In order to solve for A 2 , we rewrite Eq. (p.34) as 


= 


1 


M + m 


M± 


rn'^ aM‘^ 

M2(l — A)2 +-= (M + m)2) 


Al(0) 


vr 


(5.35) 


As indicated, in general there are at least two solutions for A 2 . In fact, because Ep also depends on 
A 2 , this equation is only implicit and actually may result in more than two solutions which have to 
be determined numerically. We have solved the implicit equation by iteration and found that below a 
critical coupling of about a = 0.4 each branch A^^^ and A 2 ^ yields exactly one solution. Above this 
coupling the solution for A^^^ ceases to exist while at the same time the branch A 2 ^ develops another 
solution. As the coupling is increased further, more solutions for A 2 ^ come into existence. 

Clearly, one needs to decide which one of these solutions corresponds to reality. For small couplings 
the decision seems relatively easy - the perturbative limit (in which A = A(0) = 1 and s = M^) of 
Eq. (1^ ) is 

, M ± m , 

■^2 = ATT— + (a) . (5.36) 

M + m 

As mentioned before, we require that in the limit of zero coupling the trial action becomes equal to 
the true action. Clearly only A^^^ in the above equation satisfies this requirement, so we shall use this 
solution up to the critical coupling of ~ 0.4. Beyond this coupling we shall contend ourselves with 
using that value of A 2 ^ which matches smoothly onto A^^^ (only one of the solutions for A 2 ^ has this 
property). 

It is not entirely clear whether this is a satisfactory choice. In particular, one is naturally led to 
speculate whether this critical coupling actually has a physical meaning - perhaps the intermediate off- 
shell propagator feels the instability of the Wick-Cutkosky model already at smaller coupling constant 
than the on-shell propagator. On the other hand, the division of solutions into A^^^ and A 2 ^ in 
Eq. ( 5.35| ) is somewhat artificial as, after all, this equation is only an implicit one for A 2 . Hence it is 
also possible that the critical coupling is just a mathematical artifact. 


5.3 Total cross section 

We now turn to the total cross section. Very close to threshold we can write 

1 


p{E,s) ~ ^A(0)( 
2 


™(+) „(-) 


1 

E 


exp 




-<3 (sthr.,a:i,Eo)C(0) 


1 


1 


.^0('Sthr.) ) .^1 ('Sthr.) ) .^0) 


(5.37) 


“Another alternative, again something which we do not pursue here, is to allow A 2 to become a function of s and 
hereby enforce all threshold positions to be precisely where they should be. 

^We remind the reader that the variational equations for the on-shell propagator cease to have real solutions for 
couplings greater than about a = 0.8, signaling the instability associated with this model. 
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where 


Tfi^ 

Eo{sthr.,xi) = mX2{M + m) - Eq 

IAq 

(5.38) 

9 

7Tl 

Ei{sthv.,xi,Eo) = mXM - Ei . 

ZAq 

(5.39) 


The square root \/E — Eq gives rise to a linear behaviour of the imaginary part of the amplitude with 
the three-momentum p in the center-of-mass system as can be seen by expanding 


s = 


^ -|- p2 + + p2 


(5.40) 


In contrast to the zeroth order variational approximation (III) where this behaviour was crucially 
dependent on the specific form of the retardation function at large values of n, in first order this now 
holds for all parameterizations. This demonstrates the increased insensitivity of the first-order result 
to the variational ansatz. At threshold one obtains 




[X2m + (A 2 - X}MY 


yJ{M + m)M X 2 m^ 


X2{M + m) - 


2 ^( 0 ) 


(am 


m ^ 


2P 


where the derivative in 


Z = \ 2 -X 2 - 


m 


aM^ d 




•p{s) 


(5.41) 


(5.42) 


m + M 11 X 2 ds 

is to be evaluated at threshold. By the optical theorem the total cross section at threshold is obtained 
from Tm A(s,0) using Eq. (|5.2D . It may readily be verified that this expression reduces to Eq. (5.7) 
in the perturbative limit. 


6. Numerical Results 

In Table 1 we collect the variational parameters determined from the variational equations at the pole 
of the propagator as well as the values for A 2 fixed at the lowest threshold. As usual we have taken 
M = 939 MeV and m = 140 MeV for the mass of the nucleon and meson, repectively. It should 
be noted that the parameters of the present ‘extended’ parameterizations are slighthly different from 
those given in Ref. (Ill) because the variation now is not constrained to give the correct (zeroth 
order) thresholds. Consequently the minimum values of the variational functional are always lower 
than previously and are nearly identical with the full variational calculation (assuming no specific 
functional form for the retardation function), except at the largest coupling constant. 

The total cross section at threshold is evaluated using the optical theorem and the analytic result 
of Eq. ( 5.41j ). The elastic cross section at threshold is determined by the real part of the scattering 
amplitude which is evaluated numerically. The ‘unitarity ratio’ crei/c’'tot is shown in Eig. ^ for the 
allowed range of couplings. For comparison we include the results obtained in first order perturbation 
theory. As can be seen the variational result remains close to the unity whereas, as is well known, 
perturbation theory strongly violates it for large couplings. 
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a = 0.1 

a = 0.2 

a = 0.3 

a = 0.4 

a = 0.5 

a = 0.6 

a = 0.7 

a = 0.8 

2/1 


1.01142 

1.01950 

1.02899 

1.04023 

1.05405 

1.07171 

1.09755 

1.15359 

2/2 


2.0702 

2.3960 

2.6302 

2.8824 

3.0850 

3.1742 

3.2194 

3.1266 

y/wi 

[MeV] 

636.09 

612.10 

587.31 

560.04 

530.05 

496.02 

453.33 

381.08 

\/W2 

[MeV] 

376.68 

370.04 

358.92 

345.96 

329.66 

308.04 

280.66 

233.72 

A 


0.97297 

0.94390 

0.91223 

0.87718 

0.83739 

0.79032 

0.72977 

0.62361 

A 2 


0.97690 

0.95215 

0.92534 

0.89586 

0.86273 

0.82414 

0.77587 

0.69838 

A(0) 


1.01508 

1.03213 

1.05177 

1.07489 

1.10312 

1.13951 

1.19203 

1.30393 

Z 


0.96087 

0.91918 

0.87429 

0.82524 

0.77045 

0.70693 

0.62748 

0.49509 

(r^)o [fm] 

0.0182 

0.0258 

0.0313 

0.0352 

0.0368 

0.0327 

-0.0164 

-0.0886 

[fm] 

0.0481 

0.0697 

0.0879 

0.1049 

0.1220 

0.1404 

0.1623 

0.1983 


Table 1 ; The variational parameters for the‘ extended’ retardation function writing Ci = g‘^yi/{327T‘^) 
and C 2 = ^Txj2 and some derived quantities. A negative zeroth order rms-radius 

indicates that < 0. 
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Figure 4; The unitarity ratio iTei/cTtot at threshold as a function of the coupling a 


g‘^j{A7rM‘^) . 


22 































We also investigated the momentum dependence of the total cross section above the elastic, but 
below the first inelastic, threshold. For very small momenta one can apply Eq. (5.2£) to obtain 
an analytic expression for the imaginary part of the forward scattering amplitude as a function of 
momentum as was done at threshold. In Fig. ^ we show the results for a = 0.8. The results for 
the three lowest momenta are calculated analytically and they join nicely to the results obtained 
numerically at larger momenta. In contrast to what was observed for the vertex function in Ref. 
(IV), here one finds a smooth behaviour of the total cross section above threshold as a function of the 
momentum. 
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Figure 5: The total cross section as a function of the center-of-mass momentum |p| at a = 0.8 . 


Conclusions and Outlook 


In this work Feynman’s variational method originally applied to the polaron problem, a non-relativistic 
field theory, was successfully generalized to study particle scattering in a relativistic setting. A central 
feature of the approach is that the variational principle is applied in the particle (or worldline) rep¬ 
resentation of field theory, making it rather similar to the non-relativistic theory. This is particularly 
so for the scalar theory under discussion in this paper - it is fair to characterise it as essentially 
Feynman’s polaron in 4 -|- 1 rather than 3 -|- I dimensions (with, in this case, phonons of non-zero 
mass). 

The method amounts to a well-dehned non-perturbative expansion of the Green functions of the 
theory. The leading term, which is the one calculated in this paper, obeys a variational principle. 


Higher order terms could be calculated as they have been for the polaron [^|, if the need arose. 
Hence, roughly speaking, the “expansion parameter” is a dynamically determined quantity which 
adjusts itself in order to minimize the corrections due to these higher order terms. It is, of course, 
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crucial, that the starting point - the trial action - contains as much of the physics of the problem as 
possible. Although practically we are limited to quadratic (but nonlocal) trial actions there seems to 
be sufficient fiexibility in the ansatz to achieve this goal. In the present work we even have generalized 
the previous trial action for better off-shell propagation without destroying the advantageous features 
of our method; the variational equations remain the same as previously, are independent of the number 
of external mesons and can be easily solved in Euclidean space. 

In any expansion scheme such as this it is desirable that the lowest order term(s) are not only 
numerically reliable, but that the expansion also exhibits some of the analytic features of the full 
theory. A trivial illustration of this statement is relativistic invariance; both perturbation theory, 
for example, as well as the method outlined here are written in a manifestly invariant form. Also, 
crossing symmetry is maintained in both approaches. Another example, one that perturbation theory 
does not share with the present approach, is that the lowest order term already exhibits all the cuts 
and poles of the exact Green function. 

A property of the theory which is connected with this, and which is maintained (at finite order) 
neither in the variational calculation nor in perturbation theory is the unitarity of the theory. In 
particular, the imaginary part of the forward scattering amplitude is connected to the total cross 
section via the optical theorem. Perturbation theory violates this very badly, something which is of 
course most obvious at tree level - here the imaginary part of the forward scattering cross section is 
identically zero because of the lack of analytic structure mentioned above. Because the variational 
calculation does not suffer from this deficiency one might suspect that unitarity is partially restored. 
The numerical results presented in this paper confirm this. 

Although this is an important result in itself, restoration of unitarity was of course not the 
prime motivation for the work presented here. Rather, we were interested in the intrinsically non- 
perturbative nature of the variational method. For this reason we concentrated our attention on a 
kinematic region where non-perturbative effects are usually most important, namely near a thresh¬ 
old. We chose the lowest theshold because this did not require a calculation of higher-point Green 
functions. 

The threshold region is important from another perspective as well. This region is sensitive to the 
infrared behaviour of the theory, while the fact that the theory is a divergent field theory means that 
ultraviolet physics plays a dominant role as well. The variational method is able to deal sensibly with 
both these regions simultaneously through the use of the prohle function A{E) in the trial action, which 
is complementary to the retardation function f{cr). Specifically, the ultraviolet (infrared) behaviour is 
simulated by the small- (large-) a behaviour of this function which corresponds to the large- (small-) 
E behaviour of A(E). 

A difficulty with investigating the behaviour of the Green function near threshold is that although 
the variational approximation to the forward scattering amplitude does indeed contain cuts which 
one can identify with one, two etc. meson production in the intermediate state, the thresholds for 
these cuts are in general only approximately equal to their physical values. The solution of this 
problem which we chose illustrates another important feature of the variational method; although 
some properties of the exact Green function are reproduced for arbitrary trial actions, others may 
be enforced by restricting the form of the trial action. In particular, in the present work the new 
parameter A 2 , associated with the internal nucleon line, was fixed so that the lowest threshold was 
identical to the physical one. 
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This is certainly a legitimate procedure, although it would be more satisfying (and one would 
expect the numerical results to improve) if the variational principle itself would provide the correct 
thresholds. This would mean, in the context of our trial action, treating A 2 as a variational parameter 
with the daunting task of solving the variational equations for this parameter. We have chosen not to 
go into this direction but rather concentrate our efforts on a generalization of the present approach 
to a more realistic field theory, namely Quantum Electrodynamics. 


Acknowledgement: We thank Nadia Fettes for a careful reading of the manuscript. 
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